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Abstract 

We consider nuclear spin-lattice relaxation rate resulted from a fractional diffu- 
sion equation for anomalous rotational wobbling in a cone. The mechanism of re- 
laxation is assumed to be due to dipole-dipole interaction of nuclear spins and 
is treated within the framework of the standard Bloemberger, Purcell, Pound - 
Solomon scheme. We consider the general case of arbitrary orientation of the cone 
axis relative the magnetic field. The BPP-Solomon scheme is shown to remain valid 
for systems with the distribution of the cone axes depending only on the tilt relative 
the magnetic field but otherwise being isotropic. We consider the case of random 
isotropic orientation of cone axes relative the magnetic field taking place in powders. 
Also we consider the case of their predominant orientation along or opposite the 
magnetic field and that of their predominant orientation transverse to the magnetic 
field which may be relevant for, e.g., liquid crystals. Besides we treat in details the 
model case of the cone axis directed along the magnetic field. The latter provides 
direct comparison of the limiting case of our formulas with the textbook formulas 
for ordinary isotropic rotational diffusion. We show that the present model enables 
one to obtain naturally the well known power law for Larmor frequency dependence 
of the spin-lattice relaxation rate. The latter is observed in some complex systems. 
From this law the dependence of the fractional diffusion coefficient on the fractional 
index is obtained to have a rather simple functional form. The dependence of the 
spin-lattice relaxation rate on the cone half- width for the case of ordinary rotational 
diffusion yields results similar to those predicted by the model-free approach. 
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1 Introduction 



Diffusometry and relaxometry is a traditional and well developed branch of 

nmr m, 0, 0, @], 0, 0, 0, 0, 0, ma, nu, ng, m, ng, ng, nu. 

The so-called fractional diffusion equation (FDE) became a new theoretical 
tool in this field pioneered by Kimmich and coauthors [13J more than a decade 
ago. Application of fractional calculus (for review see [18], [19], [20], [21], [22] 
and refs. therein) to description of fractional (or else anomalous) diffusion 
has been known for a long time (for review see [21], [25], [27], [22] and refs. 
therein). In the present paper we deal with the particular case of the so called 
sub-diffusion regime. Sub-diffusion is ubiquitous in nature and it originates in 
any fractal media due to the presence of dead ends on current ways [21] , [26] , 
[27], |22j. The most straightforward way to obtain sub-diffusion mathemati- 
cally is to generalize the ordinary diffusion equation by replacing the ordinary 
derivative in time by the fractional one of the order a. (0 < a. < 1). The result- 
ing FDE is a convenient mathematical tool although a phenomenological one 
due to the fact that there is no lucid and commonly accepted physical meaning 
of a fractional derivative at present [2TJ, [26], [27], [22] . Some theoretical mo- 
tivation for introducing the FDE arises from its intrinsic relationship with the 
so called continuous time random walk theory. In fact the FDE can be derived 
from a generalized Langevin equation with the memory kernel accounting for 
a power- law waiting time statistics of the trapping events [26] . 

Much work is carried out at present in the field of fractional Brownian dy- 
namics in proteins and liquids by Kneller, Hinsen and their coauthors [23], 
[21], [25]. By now fractional calculus finds numerous applications in NMR dif- 
fusometry and relaxometry [13], [TJ], [15], [16] . In particular it is used for the 
description of NMR relaxation in proteins |2H], [2H]- The systematic theory of 
fractional calculus in NMR is developed at present mainly due to the efforts 
of Magin and coauthors [3U], [21], [32]. However many important theoretical 
issues still remain open. For instance much studies were devoted to the case 
of translational anomalous diffusion (see [16] and refs. therein). In particular 
the generalization of the well known formula VIII. 114 from |33j for spectral 
density of the correlation function for ordinary translational diffusion to the 
case of fractional one was obtained [33], [35]. At the same time analogous 
generalization for the case of rotational diffusion has still been lacking in the 
literature. Although that for free rotational diffusion can be obtained trivially 
(basing, e.g., on the results of [36], [37] , [38], [39]) it is generally of no interest 
for applications. In practice one is encountered with some sort of restricted 
rotational diffusion, e.g., wobbling in a cone treated rigorously by Wang and 
Pecora [ID] . The latter model is widely used for interpreting NMR relaxation 
data 0, [H], [12J. The aim of the present paper is to develop the corre- 
sponding theory for fractional rotational diffusion in a cone. The approach is 
discussed mainly for practically important case of an isolated two-spin sys- 
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tern comprising a hetero-nuclear pair of non-identical spins, e.g., 15 N — H in 
protein backbone or 13 C — H in protein side chains. However the results for 
homo-nuclear spin pair are also presented for the sake of completeness. From 
the point of view of theoretical perspectives the present paper is conceived to 
attain similar goals to those of [29]. Nevertheless the implementation is basi- 
cally different. The authors of [29] consider the fractional Ornstein-Uhlenbeck 
process (sub-diffusion in a harmonic potential) and make use of the results for 
one- dimensional translational FDE to the case of rotational one. In contrast 
we consider fractional rotational diffusion in a cone within the framework of 
Wang and Pecora model that enables us to deal with the rotational FDE. In 
our opinion both models have their merits and drawbacks. On the one hand 
the harmonic potential seems to be more realistic model for restricted rotation 
than the cone with non-permeable boundary. On the other hand the authors 
of [2H] have to assume that the internuclear vectors "fluctuate only moderately 
about their average direction in a molecular- fixed frame". For wobbling in a 
cone we have no such restriction and we can extend the cone half-width 9q up 
to the limit of isotropic rotation Oq — > 7r to verify the coincidence of the results 
obtained with the known formulas. Besides the Wang and Pecora model is ex- 
actly tractable both for ordinary diffusion and for anomalous sub-diffusion. 
That is why the results of rigorous treatment of the rotational FDE can serve 
as a touchstone for approximations by necessity invoked to in the case of more 
realistic potentials. 

Experimental data that can be treated with the help of the theory developed in 
the present paper concern the magnetic field strength (or equivalently Larmor 
frequency) dependence of spin-lattice relaxation rate. In complex systems the 
latter usually exhibits a power law (1/Ti) (u) = A uj-P where < < 1 and A 
is a constant. Such dependence is observed for proteins and homopolypeptides 

0, 0, [33], [33], jsg, |5I], HZ], [22], [52], [53], 155], & [53, tissues [SB], 

nematic liquid crystals [59], [T7], liquids in porous glass [7], [T7], polymers 
[T7] , etc. Purely empirical way to obtain it is to introduce a distribution of 
the correlation times in the system, e.g., Cole-Davidson or Cole-Cole ones (see 
e.g., [42] . [45] , [60] and refs. therein). Usually this is done within the framework 
of the so-called model- free approach [30], 03, [33], [33], [3E], [32]- The generic 
feature of any such approach is that it does not touch upon the exponential 
character of the decay of the correlation function. There are more involved 
approaches to the problem that cast doubt on the validity of such simple basic 
behavior. For instance Korb and Bryant invented an ingenious and profound 
method to describe the above mentioned power law of spin-lattice relaxation 
rate in proteins [50], [55] and nematic liquid crystals [5H]- It is based on the 
relaxation mechanism stipulated by spin-phonon coupling and a postulate that 
the protons of protein structure form a percolation network with some fractal 
dimensionality df. However this model is a matter of controversy and there is 
a competing approach to analyze experimental data [61J. 
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There is profound inherent conceptual relationship between fractals and frac- 
tional relaxation [21], [26], [27], [22]. That is why it seems useful and interesting 
to explore the problem of the power law for the spin-lattice relaxation rate in 
the "language" of fractional calculus. We show that the application of the sin- 
gle order FDE enables one to obtain the desired power law for the spin-lattice 
relaxation rate (1/Ti) (cu) oc u^ 13 . Our approach leads to replacement of De- 
bye relaxation (exponential function) by fractional relaxation (Mittag-Leffler 
function), i.e., to that of Lorentz functions by those 

T , . 2rsin(W2) (ur)"- 1 

J\W) OC ; — r ; r— - 

v ' 1 + 2cos(7ra/2) (cur) a + (wr) 2a 



This spectral density leads in the extreme narrowing limit lot << 1 to the 
required power law for the spin-lattice relaxation rate with /3 = 1 — a. It is 
the well known Cole-Cole spectral density [23], [21], [25]. The fractional index 
a phenomenologically takes into account that the character of rotational dif- 
fusion may be different depending on the environment of the two-spin system, 
e.g., close packing in protein interior or loose packing at protein surface, or the 
influence of the environment depending the hydration extent of the protein, 
pressure, etc. Our approach to rotational diffusion in application to proteins 
may be regarded as a complementary one to that of [50], [55]. It enables to look 
on the problem of spin-lattice relaxation in proteins from somewhat different 
angle compared to the theory of Korb and Bryant. The approach is based on a 
simple and visual physical model of rotational diffusion in a cone rather than 
on ad hoc invented spatial distribution (fractal network) of protons in protein 
structure. Application of the model to different complex systems requires only 
the appropriate choice of the cone half-width and the value of the fractional 
index rather than profound modification of the structure of the theory. Our 
approach is in coherence with that of the paper [29J though as was mentioned 
above the implementation is basically different. 

Also the advantage of fractional dynamics approach is in the fact that it en- 
ables one to reveal links with experiments and results of computer modeling 
beyond NMR. Prior discussing them we recall (with may be somewhat un- 
conventional names) that in physics one usually distinguishes several types of 
relaxation: 1. fast relaxation (described by the exponential function), 2. slow 
relaxation (described by the famous stretched exponential or else Kolraush- 
Williams- Watts function), 3. superslow relaxation (described by a power law 
function), 4. ultraslow relaxation (described by a logarithmic function), etc. 
Remarkably that FDE enables one to obtain naturally any of the listed types 
of relaxation. The single order FDE (where the fractional index a is in the 
range < a < 1 for the case of sub-diffusion that typically takes place in con- 
densed matter problems) yields the so-called Mittag-Leffler function whose 
intermediate regime is exactly the stretched exponential function. At asymp- 
totically large time the stretched exponential behavior of the Mittag-Leffler 
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function is replaced by a power law decay. The particular case a = 1 yields 
the normal (Debye) relaxation described by the exponential function. Finally 
a distributed order FDE is known to yield the ultraslow (logarithmic) relax- 
ation at appropriate choice of the distribution p(a) over the fractional index a. 
The latter case seems to be relevant for adequate description of complex het- 
erogeneous systems such as proteins. Indeed namely the logarithmic decay is 
observed in single-particle relaxation of hydrated lysozyme powder [62]. Nev- 
ertheless for the sake of defmiteness in the present paper we restrict ourselves 
by the case of single order FDE and leave that of distributed order FDE for 
future work. Numerous interrelationships of fractional dynamics with results 
of molecular dynamics simulations are presented in [23], [23], [25], [28], [29J. 

The paper is organized as follows. In Sec. 2 the single order FDE for restricted 
rotational diffusion in a cone is used to derive the joint probability density 
function. In Sec. 3 the latter is used to obtain the spectral densities of cor- 
relation functions for dipole-dipole interaction within the framework of the 
standard Bloemberger, Purcell, Pound (BPP) - Solomon scheme. In Sec. 4 the 
dependence of the fractional diffusion coefficient on the fractional index a 
is obtained from the experimental power law for spin-lattice relaxation rate. 
Sec. 5-6 deal with the particular case that the cone axis is directed along the 
magnetic field. In Sec. 5 the spin-lattice relaxation rate for hetero-nuclear spin 
pair is obtained. In Sec. 6 that for homo-nuclear spin pair is obtained. In Sec. 7 
the general case of arbitrary orientation of the cone axis relative the magnetic 
field is considered. These results are applied to the case of isotropic random 
orientation (unweighted average) of cone axes relative the laboratory frame. 
Also a model example of predominant orientation of cone axes along or oppo- 
site the magnetic field and that of their predominant orientation transverse to 
the magnetic field which are considered. In Sec. 8 the results are discussed and 
the conclusions are summarized. In Appendix A some known mathematical 
formulas are collected for the convenience of the reader. In Appendix B some 
technical details of calculations are presented. 



2 Single order fractional case for restricted rotational diffusion in 
a cone 

We choose a laboratory fixed frame so that its z axis of the Cartesian frame 
x, y, z is directed along the constant magnetic field. The random functions F™ 
of relative positions of two spins specified below (see (JHJ)) are defined in the 
corresponding spherical frame 9, <p given by the polar angle 6 (counted from 
the z axis) and azimuthal one <j). We consider a general case that the cone axis 
is tilted at an angle ip relative the magnetic field. We direct the z' axis of the 
dashed Cartesian frame x', y', z' along the cone axis. The correlation function 
for wobbling in this cone we define in the corresponding spherical frame 6', 4>' 
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given by the polar angle 6' (counted from the z' axis) and azimuthal one 
0'. Following [ID] we consider a rod with the orientation specified by a unit 
vector u directed along its axis with spherical polar coordinates Q' = (8',<f)'). 
In accordance with |33J we following Debye assume that the rotation of the 
rod can be considered as that of the hard sphere with radius a (a is the 
length of the rod) in a medium with viscosity rj. For the ordinary diffusion in 
a cone model the rod is allowed to diffuse freely within an empty cone with 
a maximum polar angle 6' = 6q. The symmetry axis of the cone is taken to 
be the z' axis. For the diffusion in a cone model, the polar angle is restricted 
(0 < 6' < 6> ) but the azimuthal angle is not (0 < <p' < 2vr). 

Our aim is to consider sub-diffusion for rotational motion in a cone. The 
probability density for finding the rod oriented in u at time t, i. e., ^/(it,t), 
obeys the single order a FDE for rotational motion [36], [37], [38], [39] 

^)^ oA -AW(M) (1) 



where oDl~ a is the Riemann-Liouville fractional integral (see [18], [19], [20] . 
[21] . [21] . [26] . [27] and refs. therein), C a is the fractional diffusion coefficient 
(FDC) for rotation and A is the angular part of the Laplace operator in polar 
spherical coordinates 



sin^ 9' 



. a , d ( . Ql d \ d 2 



(2) 



The FDC for rotation C a has the dimension cm 2 /s a and for the case a = 1 
coincides with the Stokes formula 

c = |^ (3) 



where is the Boltzman constant and T is the temperature. 



Following the combined strategy of [ID] on the one hand and [36] 
[39] on the other hand we can write the solution of ([1]) as follows 

= 

oo oo 

E E E a \-u-(u- + l)C a \t\ a /a 2 ]Y^ * (Q'(0))Y^ (Q'(t)) (4) 



n=l m=— oo 



where E a (z) is the Mittag-Leffler function (see p], PSj, [20], [2T], [21], [26], 
[27] and refs. therein), (f2') is the generalized spherical harmonics of de- 
gree v™ [40J, the symbol * indicates the complex conjugate and the values of 
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v™ are determined by the boundary conditions on 9' defined by our diffusion 
in a cone model. The boundary condition says that there is no net change of 
the probability density at the boundary of the cone, i.e., 



89' 



(5) 



The values v™ are known functions of cos#o [10], [63]. They satisfy the re- 
quirement v~ m = u™ . The index n is defined such that v™ < v™ < v™ < ... 
for any m. The detailed calculations of u™ are presented in the tables |63j . 
The values of v™ for n > 1 increase sharply with the decrease of the confining 
volume. 

The solution PJ is subjected to the initial condition 

*(«, 0) = 5 (pf - fl'(0)) = 5 (cos 6' - cos 9'(0)) 5 (0' - 0'(O)) (6) 



The joint probability of finding the rod with orientation u(0) in solid angle 
dQ'(0) at time and u(t) in dQ'(t) at time t can be written as 



X 



2tt(1 - cos# ) 

E E E a -v™{y™ + l)C a | t | a /a 2 l ^ * (fi'(0)) ^ (7) 



oo oo 



n=l m=— oo 



The latter satisfies the normalization condition 

p(n'(t),t; n'(o),o)dn'(t)dn'(o) = 1 



(8) 



cone cone 



where the angular integrals are taken only over the volume of the cone. 



3 NMR framework for single order fractional rotation 



At the beginning of this Sec. we recall the main facts from the general the- 
ory of spin-lattice relaxation by dipole-dipole interaction suggested by BPP 
and developed by Solomon [M]. The BPP-Solomon scheme is substantiated 
by more stringent Redfield's theory (see [33] for detailed presentation). This 
scheme is developed in the frame whose z axis is directed along the constant 
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magnetic field. For the case of identical spins / the contribution to the spin- 
lattice relaxation rate constant due to rotational diffusion with the spectral 
density at a Larmor frequency Ml of the correlation function for spherical 
harmonics has the form (see VIII. 76 in [55] ) 

= 3 ™ + 1 > {J m M + J ,2) (^D} (9) 

where 7 is the gyromagnetic ratio of the nucleus, I is their spin and h is the 
Planck constant. For non-identical spins I and S we have four equations (see 



7/ 2 7^(/ + l)x 

{-^ (0) (-i--f) + ^ (2) (-i + -f)} (10) 

Here only two equations are written out explicitly because the other two can 
be obtained from them by mere changing of indexes [33]. To find the spectral 
densities J(°>(lo), J^'(uj) and J^ 2 \u) we need to know the correlation functions 
G®(t) where i = 0,1,2 . 

As was stressed in the previous Sec. we set 9' and 0' to be polar angles defining 
the direction of the axis connecting protons and ^/(9',4>',t) = \1/ (fl',t) to be 
the probability of the orientation of this axis in the direction Q' at time t. In 
the general case the axis of the cone can be tilted relative the magnetic field 
(z axis of the laboratory fixed frame) at an arbitrary angle ip. That is why 
at application to a realistic system the correlation function of internal motion 
in the cone < F^'(0)F^'(t) > internal has to be averaged over the orientations 
of cone axes with some overall distribution of the angles f(ip) characterizing 
the system of interest. That is the correlation function G^ l \t) whose spectral 
densities are to be substituted into (0) or ffTUl) has the form 

GW(t) =« F®(0)F®(t) > m ternal>overall = 
1 ^ 

-J dijj sin</> /(V) < F«(0)F«(t) >internal (11) 




VIII.88 in [23]) 

(l/Tl 1 ) = 

V ' 1 /rotat 

and 

(l/lf) = 

V ' 1 /rotat 
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To put it differently we assume that the distribution function f(if), A, u) (where 
if>,\,w are Euler angles for rotation of the dashed Cartesian frame x',y',z' 
relative the laboratory fixed one x,y,z) depends only on the Euler angle if), 
i.e., f(ip,\,u) = f(if))- It will be shown in Sec. 7 that only in this case the 
overall averaging provides the absence of cross- correlational functions with 
q ^ q', i.e., 

« F«(0)FW(t) >internal>overall= S qg >G^(t) (12) 



where 5 nm is the Kronecker symbol. The latter requirement is crucial for the 
validity of the BPP-Solomon scheme [53] • Further we consider four particular 
cases. 

a) . The cone axis is directed along the magnetic field for all cones in the system. 
It means that f(if>) = S(if>) where 5(x) is a Dirac 5— function (see Appendix 
B for technical details). This case is of little practical significance. However it 
provides direct comparison of the limiting case of our results with the textbook 
formulas from [33]. Thus it serves as a test for the validity of the present 
approach from the theoretical side. Besides the formulas obtained in this case 
without superfluous complexities are further used in more involved cases as 
building blocks. That is why we denote the correlation functions G^'(t) for 
this basic ones g®(t) = (£)/(v>)=5(V0 ■ This case is considered in 
details in Sec.5-Sec.6. 

b) . For the particular case of random isotropic distribution (unweighted aver- 
age) of cone axes relative the magnetic field we have f(if)) = 1. This case is 
considered in Sec. 7. 

c) . As an example of the case for the cone axes to be predominantly oriented 
along or opposite the magnetic field we consider the model function f(if)) = 
3 cos 2 if). This case is considered in Sec. 7. 

d) . As an example of the case for the cone axes to be predominantly ori- 
ented transverse to the magnetic field we consider the model function f(ip) = 
3/2 sin 2 ip. This case is considered in Sec. 7. 

From now and up to the end of Sec. 6 we consider the case a)., i.e., f(ip) = 
S(if>) where 5(x) is a Dirac 5— function. In this case each cone axis is directed 
along the magnetic field and we need not distinguish the dashed Cartesian 
frame x',y',z' from the laboratory fixed one x,y,z. To retain the uniformity 
of designations for correlation function of wobbling in a cone we further use 
the dashed Cartesian frame and correspondingly the dashed spherical frame 
9',(f>'. We start from the general expression for the correlation functions of 
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arbitrary order (see VIII. 13 in [33]) that in our case takes the form 

g W(t) = G^\t) m=m =< F®(0)F®(t) >internal= 

F {i) * (n'(t)) F {i) (n'(o)) P (n'(t),t; n'(o), o) dn'(t)dn'(o) (13) 



cone cone 



where i = 0,1,2 and the angular integrals are taken only over the volume of 
the cone. We need the correlation functions g^(t), g^(t) and g^ 2 \t) in order 
to substitute their spectral densities in (jUJ) or ffTOl . For our case of dipole- 
dipole interaction of two spins separated by the distance b they are defined by 
random functions F^°\ F^ and F^ [33] whose relationship with associated 
Legendre functions (cos6) is known (see, e.g., Appendix c. in [55] ) 

F (o, (fi) = ^Wf = _| pf(cos()) 

F« (Si) = *M^pvM = J_p« (cos()) exp W ) 

pm (n) = su^e*p(2^) = _l_ pf (cQ8 e) exp( . 20) (M) 



We stress once more that in our particular case f{rj}) = 8(ip) the F® (O) in the 
laboratory fixed frame are identical to F^ (Cl 1 ) in the dashed (cone-related) 
frame to be substituted in ({TBI . 

Now we have to substitute © and (|T4|) into (fTBl . We denote 

/x = cos#' (15) 

so that /io = cos #o and introduce the associated Legendre functions pjjm (/x) 



Y ^ ^ = \/^M exp^')^ (/x) (16) 



which satisfy the orthogonality properties 

' d/x PjJ (/x) P^ (/x) = flg 1 ^ (17) 



mo 



where £„ jm is the Kronecker symbol (5 n m = 1 if n = m and 5„ jm = other- 
wise). 
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Making use of 1.12.1.12 and 1.12.1.9 from |66j (see Appendix A) respectively 
we obtain after straightforward calculations 



j_+/fo 



5> a -v\{vl + l)C a | t T /a 2 



n=l 



rr(l) 
tin 



1 



K + 3) 2 K-2)^ 



{ [(^ + 3)/x 2 - l] P« M - ^/ioP^+i M}' (18) 



and 



9 (2) (t) 



(l-/i ) 2 (l + /i ) 



3 oo 



b 6 



:^+3) 2 k-2)2 



ft 



n=l 
2 



H, 



(2) 



(19) 



The calculation of G^(t) requires the formula from [?0 



K° n = Jdfi (3/i 2 - 1)P® = Az (l - 6z + Qzl) F (-z/°, v° n + 1; 2; z ) + 



Ml) 



3z (1 - 2^ ) F (-^°, »l + 1; 3; z ) + 2z 2 F u° n + 1; 4; 



20 



(20) 



where 



1-/^0 



(21) 



and F(a,b; c; x) is a hypergeometric function. Making use of ( 12"U|) we obtain 



1 °° r n 1 

J2E a \-iy° n (u n + l)C a \t\ a /a 2 ' 



(l-^o)& 6 ^ 



(0) 



{K} 2 (22) 



We denote 



r, 



(m) 



1/a 



i#(C + l)C„ 



(23) 



where m = 0,1,2. Making use of the known Fourier transform of the Mittag- 
Leffler function [23J we obtain for the basic spectral densities 
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and j^(ui) of g^(t), g^(t) and g^ 2 \t) respectively 



i (0) H 



2sin(7ra/2) 
(l-^o)& 6 



x 



E 

n=l 



(o) 



(0) 



a-1 



1 + 2cos(vra/2) 



(o) 

UJTn 



+ 



(0) 
UJTn 



2a 



r(o) rW 



(24) 



and 



T, 



(1) 



(1) 



a-1 



(1) 



+ 



UJTn 



la (i) 
tin 



X 



n=i 1 + 2 cos (7ra/2) 
+ (1 3 )t ( ^ ) _ 2) 2 { N + 3 Vo - l] ^ M - ^aP% +1 M}' (25) 



and 



2sin(7ra/2) 



E 



-(2) 



(2) 



a-1 



66 n=i 1 + 2 cos (vra/2) 
(l-/i ) 2 (l + ^o) 3 



(2) 



+ 



(2) 
UJT n 



2a 



(26) 



These spectral densities enable us to calculate any of the spin-lattice relaxation 
rates fl9|)- f|T0|) for the particular case that the cone axis is directed along the 
magnetic field. For the sake of saving room we do not write out them in 
the general form until the fractional diffusion coefficient C a is specified as an 
explicit function of the fractional index a. In the next Sec. we attain the latter 
goal on the basis of experimental data. 



4 Fractional diffusion coefficient 



As was posed in the introduction the main aim of the paper is to suggest the 
interpretation of the power law (1/Ti) (u) = A oj~" (where < < 1 and A 
is a constant) in the Larmor frequency dependence of spin-lattice relaxation 
rate that is ubiquitous in complex systems. It should be stressed that the 
fractional diffusion coefficient C a is generally an unknown function of the 
fractional index a. Ideally the function C a should be derived from some first 
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principles for each particular system under consideration. This purpose can 
hardly be achieved at present level of knowledge. That is why we can pose 
the following inverse problem: basing on the available experimental data for 
the spin-lattice relaxation rate we find the unknown function C a to satisfy the 
above mentioned power law. It will be shown below that for this purpose we 
should set 

a 2 . {ira\ 

C a = —sm — (27) 



where we denote the correlation time r and the rotational correlation time 

T~rotat 

87ia 3 r] _ Ana 3 !] r 

T = ~k^T Trotat = 3k^T = 6 ( } 



At a = 1 we obtain ([3]) as it must be. Any other choice of the function C a 
leads to strong dependence oc sin [k{1 — (3)/2] of the constant A in the power 
law on the parameter j3 that is ruled out by the experiment. 



5 Spin-lattice relaxation rate for nuclear pair with non-identical 
spins 



For the particular case that the cone axis is directed along the magnetic 
field we have for the spectral densities to be substituted in (I9l)-( fl0|) 
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(o). 



are: 



and j^(uS) are given by (!2lj) - fT26l) . It is worthy to note that if we 
identify S with, e.g., 15 N from a nuclear pair of non-identical spins 15 N — H 
then we actually need only the formula for (l/Tf 1 ^ ^ ^ (see fflQ]) ) for the anal- 
ysis of experimental data. That is why further we restrict ourselves only with 
explicit writing out the formula for this quantity. The substitution of (|27|) into 
([MD-pl and the latter into (pTOD yields 
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(29) 



This formula describes the spin-lattice relaxation rate from restricted rota- 
tional diffusion in a cone for the particular case of the cone axis to be directed 
along the magnetic field. The series in this formula is well convergent. That is 
why in practice it is sufficient to restrict oneself only by several initial terms 
in it. 

In the extreme narrowing limit bj\r << 1 the formula (1291) yields the required 
power law 



rotat 



A 



uj t l t 



(30) 



with 



= l-a 



(31) 



Here A is related to the experimental value A by obvious relationship A 
A/r 13 . The expression for A is 



{ _ Tyh 2 s h 2 S(S + l) ™ 
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(32) 



Equation (|32|) yields the relationship of the value A with the half-width of 
the cone 9 , i.e., /j, (see (|T5|) ). However as 1/™, if™ and i^ m) depend on /i 
this relationship can not be solved as an equation for the unknown # from 
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the known value of the parameter A. In practice it is only feasible to seek the 
required value of A with the help of varying # by the trial-and-error method. 

The choice of the fractional diffusion coefficient C a in the form (127)1 enables 
us to eliminate strong dependence oc sin(7ro;/2) of A on the factional index 
a. However A given by ( 1321) still exhibits some weak dependence on a. The 
latter is negligible. 



6 Spin-lattice relaxation rate for nuclear pair with identical spins 



For the case of identical spins the substitution of ff27j) into (T24|) - (|26|) and the 
latter into (Q yields 
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(33) 



It will be shown later that in the limit of ordinary (a = 1) isotropic (8q — > it) 
rotational diffusion this formula yields the well known result VIII. 105 from 



7 Arbitrary orientation of the cone axis relative the magnetic field 



In the general case of arbitrary tilted cone axis relative the magnetic field we 
need two frames (see Sec. 2 and Sec. 3). The laboratory fixed frame has the z 
axis directed along the magnetic field while the dashed (cone-related) frame 
has the z' axis directed along the cone axis. The angle between the cone axis 
and the magnetic field is ip. In ( TTBl we carry out internal averaging over the 
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rotation in the cone in the dashed (cone-related) frame. That is why we need 
the transformation of the F^' (fi) given by (Tbil) in the laboratory fixed frame 
into those in the dashed (cone-related) frame. 

As is well known a rotation of one frame relative the other is most conveniently 
described by Euler angles. We choose the angle if) as the first Euler angle 
(0 < if) < 7r). We denote two others Euler angles as A (that between the x-axis 
and the so-called iV-line (node-line) < A < 2n) and u (that between the 
iV-line and the x'-axis < u < 2tt). The formula for transformation of the 
generalized spherical harmonics at transition from the frame {0, 8} to that 
{<//, 6'} obtained by rotation of the z axis by the angle if) is [68] . [69J 

n 

P^ q) (cos 6) exp(iq<f)) = £ R { 2 s) q (iP)P^ s) (cos 9') exp [is<f)' + iqX + isu] (34) 



where 



min(2—q,2—s) 
r=max(0,—q—s) 

(2 + q)\(2-s)\ ( ip\ q+s+2r (. ^Y' q - S - 2r ( ^ 

' cos — sm — (35) 



r\(2 - q - r)\(2 - s - r)\(q + s + r)\ V 2 J V 2 



One can see that if the distribution function f(ip,X,u) characterizing a system 
of interest depends only on the angle if>, i.e., f(if>, A, u) = f(ip) then at overall 
averaging 

< ■■■ >overaii= g-j J dip J d\ J du simp f (if) , A, u) . . . (36) 





we have the factors 

J d\ exp [i(q — q')X] = 2n5 qq i (37) 
o 

J dcu exp [i(s — s')oj] = 2tt5 ss / (38) 
o 



It is namely the identity ( 13?]) that provides the applicability of the formula 
f )12p . The latter is crucial for the validity of the BPP-Solomon scheme [53] . 
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The overall averaging takes the form 
1 77 

< ••• >overall= ^ J dip Simp f(lp)... 



(39) 



and total (overall + internal) averaging is given by (jTTj) in this case. The 
function f(ip) must be normalized so that 
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(40) 
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After lengthy but straightforward calculations we obtain 



(42) 
(43) 
(44) 



where j^(co) and j^(u) are given by (j2"4"j) - (j26ft . For the case of the 

fractional diffusion coefficient (1241) consistent with the power law for the spin- 
lattice relaxation rate the explicit form of the basic spectral densities j^(co), 
j^{u)) and j^{uj) is 
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(47) 



For practical application of the theory we consider three cases. 

1. For the particular case of random isotropic distribution (unweighted aver- 
age) of cone axes relative the magnetic field we have f(ip) — 1- I n this case 
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we obtain: = 1/15; /ijjj = 1/5; /i^j = 4/5; = 4/15; /ijjj = 8/15; 
h§ = 4/5; fcgj = 2/5; Jijg = 4/15; fcgj = 16/15. 

2. As an example of the case for the cone axes to be predominantly ori- 
ented along or opposite the magnetic field we consider the model function 
f(ip) = 3 cos 2 ?/'. In this case we obtain: h$. = 1/35; h^l = 1/7; ftiS = 44/35; 

= 12/35; h$ = 28/35; h$ = 20/35; h$ = 22/35; h$ = 12/35; 

hfy = 16/35. 

3. As an example of the case for the cone axes to be predominantly ori- 
ented transverse to the magnetic field we consider the model function f(ip) = 
3/2 sin 2 ijj. In this case we obtain: ftiS = 3/35; h^l = 8/35; hS, = 20/35; 

h§ = 8/35; h§ = 14/35; h$ = 32/35; h$ = 10/35; h§ = 8/35; h$ = 
48/35. 

Making use of (|45|) - fj47|) and explicit values of h$ for these two cases enables 
us to calculate the spectral densities (j4"2l - ([4"4"I) to be inserted in BPP-Solomon 
formulas (f9~l)- (fT0l) . We do not write out explicitly the corresponding expressions 
for the spin-lattice relaxation rate to save room. However in the next Sec. we 
present the corresponding figures for spin-lattice relaxation rate for nuclear 
pair with non-identical spins for both cases. 



8 Results and discussion 

For 15 N — H nuclear pair of non-identical spins, we identify S with 15 N 
nucleus and I with H one. Thus 75 = —2712 rad s _1 Gauss' 1 and 7/ = 
26753 rad s _1 Gauss' 1 so that "is III = —0.101372. For 13 C — H nuclear pair 
75 = 6728 rad s -1 Gauss' 1 so that 75/77 = 0.251486. The righthand side in 
the formula (129]) depends on the parameters characterizing the nuclear pair 
(namely on the gyromagnetic ratios 75 and 7/ of our pair of non-identical 
spins). To plot the spin- lattice relaxation rate with the help of (I2"9"j) one has 
to choose the particular nuclear pair explicitly. That is why to be specific 
we choose the 15 iV — H nuclear pair of non-identical spins. In Fig. 1 the de- 
pendence of the spin-lattice relaxation rate on Larmor frequency in the cone 
with 6q = 55° (chosen as an illustrative example) for the case of cone axes 
directed along the magnetic field at different values of the fractional index 
a is depicted. In Fig. 2 the dependence of the spin-lattice relaxation rate on 
fractional index a in the cone with #0 — 55° for different values of Larmor 
frequency is depicted. In Fig. 3 the same for 9q = 5° is depicted as an exam- 
ple of the case of a narrow cone. In Fig. 4 the dependence of the spin-lattice 
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relaxation rate on the cone half-width 9q for the case of cone axes directed 
along the magnetic field is depicted at ult = 0.1 for different values of the 
fractional index a. In Fig. 5 the dependence of the spin-lattice relaxation rate 
on the cone half-width 8 for the case of cone axes directed along the magnetic 
field is depicted at a = 0.2 (/3 = 0.8 that is a typical experimental value) for 
different values of Larmor frequency. 

In Fig. 6 the spin-lattice relaxation rate for identical spins obtained with the 
help of (|33|) for the case of cone axes directed along the magnetic field is 
depicted as a function of the cone half-width 8q at lot = 0.1 for different 
values of the fractional index a. From this Fig. one can see that for the case of 
ordinary (a = 1) isotropic (6 = 7r) rotational diffusion in the limit of extreme 
narrowing u>lt << 1 the corresponding curve tends to the value 0.33. Thus 
the formula fl33|) yields 
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e ^ \rj 4 h 2 i(i + i) 




rotat 
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that taking into account (128]) coincides for the case with the well known for- 
mula VIII. 106 from [33] 

(1) = ?£V + 1) ££ 

\TlJ ro tat b 3k B T 



In Fig. 7 the dependence of the spin-lattice relaxation rate on the cone half- 
width 8q for the case of random isotropic distribution (unweighted average 
/(VO = 1) of cone axes relative the magnetic field is depicted at ui L r = 0.1 for 
different values of the fractional index a. In Fig. 8 the dependence of the spin- 
lattice relaxation rate on the cone half-width 6q for the case of the cone axes 
to be predominantly oriented along or opposite the magnetic field (/(VO = 
3 cos 2 ?/') is depicted at uolt = 0.1 for different values of the fractional index 
a. In Fig. 9 the dependence of the spin-lattice relaxation rate on the cone half- 
width #o for the case of the cone axes to be predominantly oriented transverse 
to the magnetic field (fijp) = 3/2sin 2 -?/>) is depicted at u>lt = 0.1 for different 
values of the fractional index a. 

As is shown in Sec. 3 in the extreme narrowing limit ojt « 1 the formula 
( 125]) yields the required behavior {\/Ti) rotat {uj) w A uj~® with — 1 — a. Such 
behavior is observed experimentally in a number of systems 03], [H], [5U] , 
[55] . [59] . [T7] (see Introduction). The typical dependence of the spin-lattice 
relaxation rate on the Larmor frequency is demonstrated by Fig.l for the 
example of the cone with the half- width 8 Q = 55°. By varying the half-width 
of the cone 9$ and the fractional index a we obtain rich behavior to fit any 
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experimentally observed power law (l/Ti) rota t(a;) ~ A uj~^ . In the extreme 
narrowing limit we get (l/Ti) rotat oc (f](T)/T) a . In the high frequency limit 
cut >> 1 we obtain from (|2"9"1) (l/Ti) rotat oc u~( 1+a ' . In this range we get 
{l/Ti) rotat oc (T/r](T)) a . The formula (I2"9"j) describes the intermediate behavior 
between these limiting cases. 

As is well known in liquids T\ ordinarily decreases with increasing viscosity, 
in some cases reaching a minimum value after which it increases with further 
increase in viscosity [33J. The variation of viscosity is caused by temperature T. 
For ordinary isotropic rotation correlation times are functions of temperature 
r n = r (r)(T)) [n(n + (see VIII. 97 in [33]) and the spin-lattice relaxation 
rate (1/Ti) ratat (t(T)) can have a maximum as a result. In the present model 
we have more options. Substitution of (|27|) into (|23|) yields 

7-H = r ( V (T)) + 1) sin (W2)r lA * (49) 

where m = 0,1,2. Thus correlation times are functions of temperature, of 
the cone half-width 9q (via the values u™) and of the fractional index a. The 
results obtained testify that variation of each of these parameters can produce 
a maximum in the corresponding dependence of the spin-lattice relaxation 
rate. For the limit case of ordinary (a = 1) isotropic (9 = tt) rotational 
diffusion at extreme narrowing condition the results obtained coincide with 
the well known textbook formula VIII. 106 from [33J. 

From Fig. 2 and Fig. 3 we see that spin-lattice relaxation rate exhibits a maxi- 
mum in the dependence of spin-lattice relaxation rate on the fractional index 
a. For a narrow cone 9$ < 30° the position of the maximum is at values of the 
fractional index a < 0.1 (see Fig. 3). As the cone becomes wider the position of 
the maximum is shifted to higher values of the fractional index a (see Fig. 2). 
For a wide cone the maximum takes place at a ~ 0.2. Our result about the 
increase of the spin-lattice relaxation rate with the decrease of the fractional 
index a is in accordance with the conclusion of the authors of [50], [55] that 
the effectively reduced dimensionality of the system (in their case it is chain 
structure of the protein compared with 3D crystalline structure) leads to more 
efficient spin-coupling relaxation. 

In Fig. 4 the dependence of the spin-lattice relaxation rate on the cone half- 
width 9 is depicted. It should be mentioned that f )29|) has uncertainties at: 
O = tt/4 because v\ = 2.0000; 9 = vr/2 because v\ = 2.0000; 9 = 3vr/4 
because v\ = 2.0000; 9 = 175° because v\ = 2.0000. However these uncer- 
tainties are isolated, and can be safely ignored. The dependence exhibits a 
maximum at #o ~ 85 4- 110°. We have a maximum at 9q ~ 85° even for the 
case a = 1, i.e., for ordinary diffusion in a cone. These maximums result from 
thorough investigation of the Wang-Pecora model [40] based on the Bauer's 
tables [S3]. The maximum for ordinary diffusion in a cone at #o ^ tt/2 is 
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similar to that given by the model-free approach. 



Indeed let us consider the most typical for practice case of the model-free 
approach when the overall motion of a macromolecule is considerably slower 
than the internal motion. In this case the expression of the model-free approach 
for the relationship of the spin-lattice relaxation rate with the order parameter 
S 2 and effective correlation time r e is given by equation (37) from jl6] 



' aS 2 + br e (l - S 2 ) (50) 



where a and b are constants independent on spatial configuration accessible for 
internuclear vector. For ordinary wobbling in a cone the relationship between 
the order parameter of the model-free approach and the cone half-width 0o is 
given by equation (A3) from 



S = I C osflo(l + cosfl ) (51) 



while that for the effective correlation time is given by equation (A4) from 



T - = 7T71 cos2 9 <> i 1 + cos i ln K 1 + cos 0o)/2] + 

D w (1 - S 2 ) { 

(1 -cos#o)/2}/[2(cos# - 1)] + (1 - cos 0„) (6 + 8 cos 9 - 

cos 2 6 - 12 cos 3 6 - 7 cos 4 # )/24 1 (52) 



where D w is the diffusion coefficient. We denote 
b 



aD. 



(53) 



Then we obtain the dependence of the spin-lattice relaxation rate on the cone 
half-width O 

— — = - cos 2 O (1 + cos0 o ) 2 + c< cos 2 O (1 + cos # ) 2 {In [(1 + cos 0„)/2] + 
a T\ 4 [ 

(1 - cos0 o )/2}/[2(cos0 o - 1)] + (1 - cos0 o )(6 + 8cos# - 

cos 2 O - 12 cos 3 o -7 cos 4 # )/24 1 (54) 
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The dependence of this spin-lattice relaxation rate on the cone half-width 6q 
is depicted in Fig. 11 at several values of the parameter c. One can see close 
qualitative similarity of the curves with our result for the case of ordinary 
wobbling in a cone. 

The BPP-Solomon scheme requires essential extent the isotropy (but not to- 
tal isotropy) for its validity. In Sec. 7 we show that in application to rotational 
diffusion in a cone it remains valid for systems with a distribution of cone 
axes depending only on the tilt relative the magnetic field but otherwise being 
isotropic. This residual isotropy provides the requirement ( )37|) that is neces- 
sary for the absence of cross-correlational functions (fl2|) . The latter in turn 
is crucial for the validity of of the BPP-Solomon scheme. We consider the 
aforesaid a so important issue that would like to reiterate it in other words 
with complete definiteness. We develop the theory for the general case of arbi- 
trary orientation of the cone axis relative the magnetic field (laboratory fixed 
frame). We show that when the cone axis is tilted at an arbitrary angle tp 
to the magnetic field but otherwise is oriented isotropically then at overall 
averaging the crucial requirement for the validity of the Bloemberger, Pur- 
cell, Pound - Solomon scheme ( I12p (that of the absence of cross-correlational 
functions with q ^ q') is retained. Thus we explicitly prove the consistency 
of combination of the textbook formulas for the BPP-Solomon scheme with 
rotational motion in a cone. 

Practical applications of our results depend on the choice of a distribution 
function for overall averaging over all orientations of the cone axis with respect 
to the laboratory fixed frame (over the angle tp). This distribution function 
f{ip) is a characteristic of the system of interest. Under the assumption of 
isotropic random orientation of cone axes relative the laboratory fixed frame 
the results obtained can be applied to powders. We consider this case of un- 
weighted average over all orientations of the cone axis with respect to the 
laboratory fixed frame (f(ip) = 1) in Sec. 7 and provide corresponding data 
for the spin-lattice relaxation rate in Fig. 7. Also in Sec. 7 we consider a model 
example of predominant orientation of the cone axis along or opposite the 
magnetic field {f{ip) = 3 cos 2 ip) and that of their predominant orientation 
transverse to the magnetic field (f{ip) = 3/2 sin 2 ip). The results for these 
cases may be relevant for, e.g., liquid crystals. For both of these cases we also 
provide corresponding data for the spin-lattice relaxation rate in Fig. 8. and 
Fig. 9 respectively. The particular case of the cone axis directed along the 
magnetic field (f(tp) = where 5(x) is a Dirac delta-function) is of lit- 
tle practical significance. However this case provides direct comparison of the 
limiting case of our formulas with the textbook formulas from [33] . Thus it 
serves as a test for the validity of our approach from the theoretical side. In 
the limit of isotropic motion (i.e., the cone half- width = it) our formulas yield 
those of [33J as it must be. 
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Fig. 7, Fig. 8 and Fig.9 show that the results for practically relevant cases of 
specific distributions of cone axes relative the magnetic field are qualitatively 
similar to those for the model case of cone axes directed along the magnetic 
field (Fig.4). However quantitatively the values of spin-lattice relaxation rates 
differ from those for the model case. For very wide cones (Oq > 100°) the 
results for three cases are rather close to each other because when the cone 
half-width is large its tilt relative the magnetic field becomes of minor impor- 
tance. Appreciable difference between three cases takes place in the region of 
moderate cone half-widths 25° < 6q < 90°. The comparison of Fig.4, Fig. 7, 
Fig. 8 and Fig.9 shows that 



Thus the more is the contribution of the orientations of the cone axes in the 
system transverse to the magnetic field the more efficiently spin-lattice relax- 
ation proceeds. To prove it explicitly we plot in Fig. 10 the the dependence of 
the spin-lattice relaxation rate on the tilt angle of the cone axis relative the 
magnetic field for the cone with the half-width 6 = 55° as an example. This 
case corresponds to the distribution function f(ip) = 5(ip — ipo) [cos (?/'o/2)]~ 1 
(see Appendix B for technical details). One can see that when the cone axis 
is transverse the magnetic field the spin-lattice relaxation proceeds most effi- 
ciently 

In closing we touch upon the most difficult and troublesome problem of the 
physical interpretation of the results obtained in the present paper. Regret- 
fully the latter is hampered by the general flaw of the fractional calculus and 
the FDE that there is no lucid and commonly accepted physical meaning of 
a fractional derivative at present. As a consequence there is no clear physical 
interpretation of the fractional index a and the results based on its varia- 
tion. Despite the aforesaid the application of the FDE to the description of 
relaxation and diffusion in complex systems is certainly a mainstream among 
modern approaches to this problem (see, e.g., monographs [21], [22] and re- 
views [22], [27J and refs. therein). Slowly this approach is finding its way into 
NMR P3], [30], [SI], [S3, [33], [3S], EE], |2S]. Notably the paper [29] is aimed 
to apply this approach to the description of rotational motion. The present 
paper is also aimed to attain it. In fact the FDE is a phenomeno logical tool 
that can find at least three different physical interpretations. 

a. The time- dependent part of its solution is a Mittag-Leffler function that 
is a natural generalization of an exponent whose spectral density is the Cole- 
Cole one [28], [29]. Thus the Mittag-Leffler function is merely a concise and 
elegant way to write out the integral over the exponents with the Cole- Cole 
distribution. The physical meaning of the fractional index in the FDE appears 
to be that of an empirical parameter of the Cole- Cole distribution that a 
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practitioner in NMR or dielectric relaxometry chooses to fit experimental data. 
The authors of [29J incline to choose this interpretation. However there are 
other options that can be used even when there is no distribution of relaxation 
times but the rod of an isolated spin pair moves on a fractal or on a lattice 
whose knots are traps with some distribution of waiting times. 

b. The FDE is well suited to describe diffusion on a fractal [26]. The fractional 
index in it is directly related to the fractal dimensionality of the system of 
interest. 

c. The FDE can be derived from the continuous time random walk theory 
of the diffusion on the lattice with traps [26J. The physical meaning of the 
fractional index is that of a parameter in the power-law distribution of the 
trapping events. 

At such abundance of possibilities to choose one of them at the level of general 
theoretical consideration (that the present paper deals with) and to insist on it 
seems to be somewhat narrow-mindedly. Of course at application to particular 
systems one can find convincing reasons to choose a specific interpretation. In 
our opinion the most honest attitude at the level of the present paper is to 
admit that there is no commonly accepted physical meaning of the fractional 
index and the results based on its variation (or more exactly there are too many 
of them). The latter statement does not mean that the FDE is inapplicable. It 
should merely be perceived properly as a phenomenological tool. This is the 
point of view accepted in the present paper. 

We refrain ourselves from further discussing this issue. The derivation of (PQ) 
is strictly speaking beyond the scope of the present paper. We feel that fur- 
ther attempts to motivate ([T]) from our side would be mere reiteration of the 
arguments presented in [23], [21], [25], [2S], [2H] within a more general con- 
text. In any case the FDE is so well suited for phenomenological description 
of relaxation processes obeying fractional power laws [21], [26], [21], [22] that 
it seems reasonable and timely to apply it to nuclear spin-lattice relaxation 
from restricted rotational diffusion and to develop a formal theory. The paper 
[2U] is a step in this direction. The present study is also a development along 
this line and in our opinion for the particular case of fractional wobbling in a 
cone the aim has been attained. It should be stressed that the absence of the 
power-law is merely a particular case of the present model (the value a = 1 
for a parameter a that can take values from the range < a < 1). Thus the 
absence of the power-law in an experiment does not invalidate the model. In 
our opinion it is useful to have a formalism that can describe well established 
results as a particular case and be easily generalized (by simple varying a sin- 
gle parameter rather than by introducing a distribution of correlation times 
via some integral) to describe more complex behavior at will. 
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We conclude that fractional wobbling in a cone is a relatively simple, ex- 
actly tractable (within the range of validity of BPP-Solomon scheme) and 
convenient model for phenomenological description of experimental systems 
exhibiting power law in Larmor frequency dependence of nuclear spin-lattice 
relaxation rate. 
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Here we present two known mathematical formulas 1.12.1.12 and 1.12.1.9 for 
the the associated Legendre function P{f (x) from the table of integrals 
The formula 1.12.1.12 is 



(x) 



:i - x 2 )^ 2 



[ dx x(l - x 2 )^ /2 P^ v , , 
{ [{y ± n + 2)x 2 - l] (x) + (jx-v- l)xP^ x (x) } 



The formula 1.12.1.9 is 

, (i _ r 2\( M +l)/2 

dx (I — x 2 Y /2 P M (x) = — — P M+1 (x) 
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The polar angle ip at operations with the distribution function f{ip) imposes 
some peculiarities in treating the case of Dirac 5— function 

cos (-00/2) 

We stress that the case of the cone angle oriented along the magnetic field 
f{ip) = considered in Sec. 4- Sec. 6 is a particular case of this distribution 
function corresponding to the value ipo = 0. First let us prove the normaliza- 
tion requirement ( HOT) . We have (taking into account that 5(2z) = S(z)/2) 

dip sin^ 6{ip-ipo) =2 J d ( — J sin - cos ^5 \2 - ^ 

^ ' ^ 
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77/2 ( v\ 1 2 
dx sin x cos x 5 \ x = 

o v 7 



J d(cosx) cosx 5 ^x — yj 



J dy y 5 ^arccosy - y^J = ^ dy y 5 \y - cos y^J 



cos — 
2 



This calculation serves as a model for operations at calculating the values of 
h^A in (j4Tj) with the distribution function = S(ip — ip ) [cos (ip /2)]~ 1 . 

The general rule takes the form 



1 



cos (i> /2) 



dif) smip 5{ip - i>o)q(if>) = 2q(ip ) 



o 



where q(ip) is an arbitrary function of if). 

At ipo = we obtain hS. = 2, fti^ = 2 and /iSj = 2 while /ij^ = at 
m n. Substitution of these values into (j42l)-(l4ll) yields J (2) (w) = j (2) (w), 
J(i)( w ) = and J(°)(a;) = j {0) (co) as it must be. 
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Fig. 1. Spin-lattice relaxation rate for 15 N — H nuclear pair of non-identical spins 
from fractional rotational diffusion for the case of cone axes directed along the 
magnetic field (eq. (|29p ) in the cone with 8q = 55° as the function of Larmor 
frequency ujl at different values of the fractional index a: a = 1 (thick line); a = 0.9; 
a = 0.8; a = 0.7 (thin line). 
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Fig. 2. Spin-lattice relaxation rate for N — H nuclear pair of non- identical spins 
from fractional rotational diffusion for the case of cone axes directed along the mag- 
netic field (eq. (|29|) ) in the cone with Oq = 55° as the function of the fractional index 
a at different values of Larmor frequency ul- ujlt = 10 -4 (thick line); oj^t = 10~ 35 ; 
UlT = l(r 3 ; oj l t = 1(T 2 - 5 ; uj l t = 1(T 2 ; oj l t = KT 1 ' 5 ; oj l t = 1CT 1 (thin line). 
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Fig. 3. Spin-lattice relaxation rate for 15 N — H nuclear pair of non-identical spins 
from fractional rotational diffusion for the case of cone axes directed along the mag- 
netic field (eq. (|29p ) in the cone with 9q = 5° as the function of the fractional index 
a at different values of Larmor frequency u>l: uj^t = 10~ 4 (thick line); ojlt = 10 -3 ' 5 ; 



uj l t = 1(T 3 ; oj l t = 1(T 2 - 5 ; uj l t = lO"*; uj l t = 10 



1-2. 



1-1.5. 



colt = 10 1 (thin line). 
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Fig. 4. Spin-lattice relaxation rate for 15 N — H nuclear pair of non- identical spins 
from fractional rotational diffusion for the case of cone axes directed along the 
magnetic field (eq. (|29|) ) as the function of the cone half-width 9q (in degrees) at 
colt = 0.1 for different values of the fractional index a: a = 1 (thin dots); a = 0.8 
; a = 0.6 ; a = 0.4 ; a = 0.2 (thick dots). 
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Fig. 5. Spin-lattice relaxation rate for 15 7V — H nuclear pair of non-identical spins 
from fractional rotational diffusion for the case of cone axes directed along the 
magnetic field (eq. (|29p ) as the function of the cone half- width 9q (in degrees) 
at a = 0.2 for different values of the Larmor frequency: w^r = 1 (thin dots); 
UlT = HT 1 ; oj l t = 1(T 2 ; uj l t = 1(T 3 ; uj l t = (thick dots). 
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Fig. 6. Spin-lattice relaxation rate for nuclear pair of identical spins from fractional 
rotational diffusion for the case of cone axes directed along the magnetic field (eq. 
(f33|) ) as the function of the cone half-width 6q (in degrees) at wlt = 0.1 for different 
values of the fractional index a: a = 1 (thin dots); a = 0.8 ; a = 0.6 ; a = 0.4 ; 
a = 0.2 (thick dots). 
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Fig. 7. Spin-lattice relaxation rate for 15 N — H nuclear pair of non-identical spins 
from fractional rotational diffusion for the case of random isotropic distribution 
(unweighted average f(ip) = 1 taking place in powders) of cone axes relative the 
magnetic field as the function of the cone half-width 9q (in degrees) at ult = 0.1 
for different values of the fractional index a: a = 1 (thin dots); a = 0.8 ; a = 0.6 ; 
a = 0.4 ; a = 0.2 (thick dots). 
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Fig. 8. Spin-lattice relaxation rate for N — H nuclear pair of non- identical spins 
from fractional rotational diffusion for the case of the cone axes to be predominantly 
oriented along or opposite the magnetic field (f(ip) = 3 cos 2 if) that may be relevant 
for, e.g., liquid crystals) as the function of the cone half- width 9q (i n degrees) at 
lolt = 0.1 for different values of the fractional index a: a = 1 (thin dots); a = 0.8 
5 ol = 0.6 j ol = 0.4 5 ol = 0.2 (thick dots). 
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Fig. 9. Spin-lattice relaxation rate for 15 N — H nuclear pair of non-identical spins 
from fractional rotational diffusion for the case of the cone axes to be predominantly 
oriented transverse to the mag netic field (f(ip) = 3/2 sin 2 V> that may be relevant 
for, e.g., liquid crystals) as the function of the cone half- width 9q (in degrees) at 
ojlt = 0.1 for different values of the fractional index a: a = 1 (thin dots); a = 0.8 
; a = 0.6 ; a = 0.4 ; a = 0.2 (thick dots). 
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Fig. 10. Spin-lattice relaxation rate for 15 iV — H nuclear pair of non-identical spins 
from fractional rotational diffusion in the cone with the half-width 9q = 55° as 
the function of the tilt angle t/'o of the cone axis relative the magnetic field (corre- 
sponding to the distribution function /(VO = —ipo) [cos (ipo/2)]" 1 ) at w^r = 0.1 
for different values of the fractional index a: a = 1 (thin line); a = 0.8; a = 0.6; 
a = 0.4; a = 0.2 (thick line). 



42 






c=50 




c=70 




c=90 




c=110 



So (degrees) 



Fig. 11. Spin-lattice relaxation rate for the model- free approach (eq. (|54j) ) as the 
function of the cone half- width 9q f° r different values of the dimensionless parameter 
c (eq. (J53J)): c = 50 (thick line); c = 70 ; c = 90 ; c = 110 (thin line). 
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